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Abstract 

We develop a general theory of a boson decomposition for both local and non-local interactions in lattice fermion models which 
allows us to describe fermionic degrees of freedom and collective charge and spin excitations on equal footing. An efficient per- 
turbation theory in the interaction of the fermionic and the bosonic degrees of freedom is constructed in so-called dual variables 
in the path-integral formalism. This theory takes into account all local correlations of fermions and collective bosonic modes and 
interpolates between itinerant and localized regimes of electrons in solids. The zero-order approximation of this theory corresponds 
to extended dynamical mean-field theory (EDMFT), a regular way to calculate nonlocal corrections to EDMFT is provided. It is 
shown that dual ladder summation gives a conserving approximation beyond EDMFT. The method is especially suitable for consid- 
eration of collective magnetic and charge excitations and allows to calculate their renormalization with respect to "bare" RPA-like 
characteristics. General expression for the plasmonic dispersion in correlated media is obtained. As an illustration it is shown that 
effective superexchange interactions in the half-filled Hubbard model can be derived within the dual-ladder approximation. 
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1. Introduction 

A self-consistent description of electrons, magnons, phonons and other fermionic and bosonic degrees of freedom 
in solids is still a challenging problem for strongly coiTelated systems where the simplest approximations are usually 
not sufficient. One of the classical examples is the theory of itinerant electron magnetism, where the d-electrons 
form both the coiTelated energy bands as well as the quasi-localized magnetic moments and the magnon spectrum 
HI. In this problem the conventional random phase approximation (RPA) gives a wrong description of the finite 
temperature effects whereas the proper corrections to RPA can be successfully calculated only in a very close vicinity 
of the Stoner instability point LL^. The main difficulty is the coexistence of localized (atomic-Uke) and itinerant 
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features in the behavior of d-electrons in solids |I3]|4]|5|. Recent progress in the theory of correlated electron systems 
is related to the development of Dynamical Mean Field Theory (DMFT) (6], which maps the problem of fermions with 
local interactions on a crystal lattice onto a self-consistent effective impurity problem. The latter can be solved very 
efficiently within Continuous Time Quantum Monte Carlo (CT-QMC) Q or exact diagonalization (ED)[6) schemes. 
As a result, both atomic-like features such as multiplet formation and band dispersion can be taken into account 
simultaneously LSj^. The combination of the DMFT approach with realistic electronic structure calculations IfTOlfTTl . 
in particular, for transition metal systems |l5]|T2l opens a new way for the description of electronic degrees of freedom 
in solids. A generalization of the DMFT approach to non-local interactions and consistent description of collective 
bosonic excitations |,10l[T3lfT4iri5 | boost first realistic applications of a so-called GW-i-DMFT approach [161 . 

Here we construct a theory, that is capable to describe systems with strongly correlated electrons interacting with 
collective bosonic modes. These modes can appear explicitly in an original Hamiltonian (for instance, phonons) or 
arise from electron-electron interactions (like plasmons or magnons). We assume that the Bose excitations cannot 
be considered as an ideal Bose gas and boson-boson correlations are important. Among numerous examples of such 
systems are the cases of weak itinerant ferromagnets 12 and, more generally, any materials near quantum phase 
transition points [17J . 

The Feynmann diagrammatic technique became the most popular tool in the quantum many body theory as a 
compact and elegant representation of a perturbation expansion in powers of interaction (781 1191 . The main advantage 
of the diagrammatic approach is an opportunity to sum up infinite series of relevant diagrams going beyond the formal 
perturbation expansion. The mainstream paradigm is that summing up more sophisticated series expands the area 
of applicability. In some cases leading series of diagrams can be separated and summed up with the use of small 
parameters which are not directly related to the interaction strength. Important examples are the Gell-Mann Bruekner 
theory of the high-density electron gas Il20l . the gaseous approximations of low density Fermi 11211 and Bose Il22l 
systems, the s-expansion in theory of critical phenomena |23|, and the expansion in the theory of magnetic 
impurities [|24I . Historically, the DMFT was introduced as a leading order approximation in a formal parameter 
l/d, where d is the dimensionality of the space. DMFT is exact in the limit of — > oo, at the same time attempts 
to construct a regular expansion in l/d turn out to be not practical f25l. Nevertheless, DMFT appears to be quite 
a good approximation for many three- and even two-dimensional systems 1.6. .10. .12,1 , since it provides reasonable 
interpolations between non-interacting and atomic limits for any dimensions. In terms of the diagrammatic approach 
the DMFT corresponds to a summation of all local in space contributions to the self-energy f6l. 

Recently we developed the dual fermion scheme ll26ll27ll which introduces new variables in a path integral gener- 
ating an alternative diagram technique such that the bare propagators (free dual fermions) correspond to DMFT. This 
means that the main local physics is taken into account by the change of variables while the new diagrams describe 
only non-local corrections to the self-energy. Interaction terms in the dual variables look more complicated than the 
original ones: the interactions become retarded and contain not only pair terms but also multi-fermion contributions. 
However, an advantage of dual fermion transformation is a faster convergence of the diagram series. For example, a 
simple ladder summation gives quite accurate description for the single band Hubbard model in strongly correlated 
regime |28|. Thus, the dual fermion approach is an example of an alternative strategy in the quantum many -body 
theory, when a proper choice of variables allows one to consider a simpler series of diagrams. Another example of 
such strategy is the use of auxiliary fermions and bosons in the Hubbard model which allows one to describe the subtle 
effect of the growth of the effective mass near metal-insulator transitions, already on the level of the mean-field ll29ll . 

The consideration of collective (bosonic) degrees of freedom in systems of interacting fermions within the con- 
ventional technique necessarily involves summations of sophisticated diagrams. The dispersion law of bosonic exci- 
tations (e.g., magnons or plasmons) is determined by the poles of the corresponding susceptibilities (the two particle 
fermionic Green's functions) and the minimal approximations which give rise to these poles, even in the case of a weak 
interaction corresponding to the ladder summation (random phase approximation, RPA [18 30 3U). In order to take 
into account magnon-magnon interactions one needs to consider more complicated, parquet-hke diagrams ^ |32l . 
For strongly correlated systems one can consider the same procedure within the dual fermion approach. Up to now 
this is only done at the level of ladder summations |28 1. As we pointed out, the above considerations of complicated 
diagrams can be replaced by a proper choice of variables in the path integral. Following this strategy we propose an 
extension of dual fermion scheme explicitly adding bosonic variables. It is important to note that the present theory 
is not a direct extension of the dual-fermion approach [26.1 . We will see that, even for the Hubbard-like models with 
on-site interactions only, it gives a different result because of the account of collective modes (magnons). 
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A simultaneous description of fermionic and bosonic correlations on the equal footing was started within the Ex- 
tended Dynamical Mean Field (EDMFT) framework [TJ, T4|. The strong-coupling generalization of EDMFT f 33l 
gives a useful extension of the Hubbard-I like perturbation expansion for lattice models [34 1. EDMFT properly in- 
cludes the physics of local fermionic and bosonic correlations. At the same time, a dispersion of collective excitations 
within EDMFT is determined by the non-local part of the interactions only. In particular, to describe the dispersive 
magnons within EDMFT one needs to explicitly add direct Heisenberg exchange interactions to the Hamiltonian. Ob- 
viously, this is not sufficient since the indirect exchange interactions (superexchange, double-exchange, RKKY, etc.) 
contribute as well to the the magnon dispersion. Moreover, typically, they give the main contribution. 

In the present paper we construct a formally exact diagram technique where the bare fermionic and bosonic 
propagators correspond to EDMFT, whereas the indirect exchange is given by diagrammatic corrections. Therefore, 
the direct and indirect interactions enter the magnon dispersion law on equal footing. 



2. Definitions 

We consider a lattice fermionic problem with the intersite hopping and generic non-local interactions. For sim- 
plicity, we present only formulae for the single-band periodic lattice, but a generalization of all expressions to a 
multiorbital case is straightforward. The action is 

S ^'^Sat[cl.,c,]+ ^ ercI^^c.+r^o- + ^ VRnP*nPr+Rn (1) 

Here c^, c are Grassmann variables, the indices r,R run over the lattice sites and translations, k stands for the wave- 
vector, T is the imaginary time, co - 2n(n + l/2)/y6 and Q = Inn/jS are, respectively, the fermionic and bosonic 
Matsubara frequencies (n is integer and /5 is the inverse temperature), and cr is the spin index. Prefactors y6 ' and ' 
are implied for sums over frequencies and wavevectors, respectively (A^ is the number of lattice sites). The "Bosonic 
variable" p in general is a vector quantity, whose components for single band models can be the charge-density or the 
spin-density operators: pj = Yj<to-' ^a-si-g-.Ca-' - p, where s-i with j = {0,x,y,z} are Pauli matrices (1, cTj., cr,,, cr^,), and 
the quantity p is introduced to satisfy the condition < p >= 0. 

To simplify the notations, we skip the index j throughout the paper, but one should bear in mind that relevant com- 
ponents of p are to be accounted in a practical calculation. For example, if the physics is determined by charge fluctu- 
ations only (as may be the case for the Coulomb interaction or an electron-electron attraction mediated by phonons). 



The Heisenberg isotropic 



only the density-density term y^^PoPo is important in the interaction part vj^^ of Eq.|Ti, 
exchange interaction corresponds to V " = V- - V"' = J, etc. An important case of theHubbard model corresponds 
to y = 0, where the spin degrees of freedom Px,y,z do not show up in the bare expression part of Eq.([T]i, but will appear 
in the following formulation, since the low-energy behaviour of system is governed by the spin fluctuations. 
Within the path integral representation the atomic action Sar is given by: 



■ ^(/w + yU)4,^C^^ + f/c|c|c|c|(iT, 



(2) 



where U is the local Coulomb interaction. 

We are interested in the single-electron Green's function G^t = - < '^/•tC^Lq^^o > and the bosonic correlator 
- - < PrrPl^i)^^!) >■ It is an important point that an account of the bosonic quantity is fundamental: we suppose 

that collective modes are relevant for the physics of the system. 



3. Local approximations 

The general idea of effective-medium approximations is to reduce the lattice problem ([T]i to the so-called impurity 
model, describing a single atom (or small cluster) in a Gaussian environment. For a system with bosonic and fermionic 
degrees of freedom, a single-site impurity problem is described by the action 

+ ^ Anp*apn- (3) 

o) n 
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Obtaining a numerical or analytical solution for the impurity problem is supposed to be much simpler We later 
assume that its Green's functions: gr — - < Ctc|^q >imp and;^f^ = - < PtP^^q >mp are known, as well as higher-order 
correlators. In practice, they can be calculated by a CT-QMC scheme fl\ or by exact diagonalization if the dimension 
of the Hilbert space is not too large. Normally x is called the full polarizability. We will call it the bosonic Green's 
function, to emphasize the formal analogy between the equations for g andx in our formalism. 

For the so-called local approximations, only the Green's functions of the impurity model gai,XQ are taken into 
account. In a nutshell, the initial problem ^ is replaced by an auxiliary lattice with Gaussian nodes; those nodes are 
characterized by the Green's functions gaj,Xn- Such a consideration gives us the well-known expressions for lattice 
Green's functions [14], which are denoted with calligraphic letters in this paper: 



= —r—. r (4) 



X, 



1 

gj + A„ - St 

1 



where k is the wave vector and Fourier transforms of the hopping and interaction parameters are introduced in the 
obvious way. One can see that the fermionic and bosonic self-energies of these expressions are A:-independent, that is, 
local in space. 

The hybridization functions A and A should be defined in such a way, that the impurity problem provides us the 
on-site (local) properties of the solution for the lattice problem. Such a requirement is physically quite reasonable. 
Also, it is closely related to Feynman's variational approach in dual fermion variables |27 |. 

Most common are so-called static and dynamic mean-field theories. In the static approach, the hybridization 
functions A and A are chosen to be frequency independent. Their values are defined from the equalities for the static 
parts of Green's functions: 

Z^""-Z^"'' 

(li (Ilk 

Z MF _\' y^ 



,MF 

OK 



In the dynamical EDMFT approach A^^ and An are adjusted so that local parts of the lattice Green's functions 
coincide with the corresponding quantities of the impurity problem/or all frequency (or time) arguments: 



Yj&cok, (6) 



K 

It is worthwhile to note that the standard definition of the EDMFT |T4l introduces an auxiliary bosonic field called 
(p, and the second line of Eq.(|6]l is re-expressed as a condition for the Green's function of that bosonic field. At the 
EDMFT level, the conditions for X and X^ are equivalent. We discuss this issue in detail in Appendix and show that 
similar (although not the same) bosonic field can be introduced in non-local approximations beyond the EDMFT. 

Some "mixed" local approximations can also make sense. For instance, if the hopping is absent (or negligible), it 
is natural to consider a scheme with A = and to adjust bosonic hybridization only. Contrary, A = corresponds to 
the "purely fermionic" DMFT scheme 0|. 

4. Dual transformation 

It is easy to make a connection between the lattice and impurity actions: 

{Sk - Aa>cr) '^Ika^ojka + ^ (VilK - ^q) PnKPnK (7) 
r bj,k,o- n,K 
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Further on, we use a general matrix form of the Hubbard-Stratonovich transformation to decouple the term 
clji^^Eaika-CojkcT with Ei^kcT - - sii into the pair of dual fermionic Grassmann-fields f^i^^,faikcr and the term 
PciKj^siKP^Kj' with W^f, = - V^^ into the pair of dual bosonic complex-fields ?7qa:;' ^^kj (a matrix notation 
is used to shorten the expressions): 

/ D [c"^ c] e^*^" = / D [c"'", c] det [a-/Ea-/) J D [/"'",/] 

(8) 

jD\p*,p] eP'^P = jD\p\p] det (ffiW-'o-i) / D [t]\ 77] e-';>/.w-'«*';+';-«*p+p>/.';, 

where aj, o-/, are arbitrary, but local in space scaling matrices. 

Note that whereas the decoupling of Grassmanian variables is mathematically straightforward, more comments 
are needed on the second line of Eq.ljSj). The first comment is that p is Hermitian, so that P*q_j^j = PnKj, and 
the integral D[p*,p] is taken in fact over the variables from a half-space. The dual variables inherit this property, 
VLa-Kj ~ ^nKj- The second comment is that the matrix ahW^^ai, is not positively defined, and integrals over D[ri*, 77] 
should be re-defined to ensure convergence. Here we describe the corresponding procedure, following, in particular, 
Refs. ll35lfT4l. 

For illustration purposes, we suppose a to be just a real constant factor; the generalization to a complex matrix 
a is straightforward. Let us introduce a basis where W"' is diagonal. Let us number the states in this basis with an 
index /, and consider some matrix element W//. For a positive real part of Wu, the corresponding part of the action can 
be decoupled as follows: 

where we explicitly introduced the Hermitian and anti-Hermitian parts of p = p' -1- ip". One can see that for 77J and tj" 
being interpreted as real and imaginary parts of a complex number 77, this decoupling becomes an explicit form of the 
second line of Eq. (j8]l, written for a particular state I. The integral limits for positive Wu can be chosen along the real 
axis, from -00 to +00 for both 77J and 77". 

The key observation is that for Wn having a negative real part, the decoupling (j9| still holds, up to the integration 
path along imaginary axis, from -ioo to +100 for 77J and in the opposite direction, that is +ico to -ico, for 77". Thus, 
the integration over D[ri*,ri] can be considered as a symbolic notation for the integration over D[?7',77"] with 77 - 
77' + irj", rf - T]' — irj" substituted in the integrand. The integration path over each component of vectors 77', 77" 
is explicitly defined for the basis where aWa is diagonal: the integration goes along the real or imaginary axis, 
corresponding to the sign of the real part of the diagonal matrix element Wn. 

One can check that such a definition of J D[ri* , 77] still allows one to use the standard algebraic manipulations of 
path integrals. In particular, to obtain the Green's function of dual bosonic variables (5 = -/ < 77 ■ 77* > one uses the 
standsrd procedure of differentiation with respect to o-W 'a, as can be explicitly shown from Eq. On the other 
hand, the statistical properties of dual variables can be unusual because of the imaginary axis integration path. For 
instance, similarly to the fermionic dual Green's function EtII . < 77 ■ 77* >=< rj' ■ rj' > + < rj" ■ rj" > is not necessary 
positive-defined and consequently the imaginary part of its Fourier transform is not always negative. 

Let us continue with the derivation of the dual formalism. The Hubbard-Stratonovich decoupling in the partition 
function with the action ([T]) gives (in matrix form) 

/ D[c^'c]e-^rSi,„AcUcr]+cHA-s)c+p'(A-v)p ^ det(Q'-i(A - £)a^') det(a'fe(A - V)-^ab)x 

(10) 

X J D[c^c] J / D[77*?7]e"^''^""''^'^"'^'^^-'^'"'^'^^'^'"^-^^''*"*''^''*'^''''"-^^ 
Integrating out the "localized" fermionic Grassmann variables c' , c with a = g^^ for the fermionic and a = Xn for 
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(11) 



the bosonic part, we arrive at the following action in dual variables 

'^nK = (Xn + An - VkT^ -xn 

where y is a full four-point fermionic vertex of the impurity model 

r^-it-i'n = (12) 

goj+ngio'-ngojgcu' 

{6 is the Kronecker symbol), and /i is a "mixed" quantity 

- < c^+aclaPn >imp - < P >imp guSn 

'iojn - ■ (13) 

gcjgoj+nXn 

The second term in the nominator typically equals zero as < p >„„p vanishes. The values of y and A are related, since 

_ V t 

Pn - luoj' c^'o-Vo-'C^j'-ao-s 



1- X y''>^<^'oC''^<T<T'fu,'-ii\ (14) 

It is important to note that for the Gaussian impurity problem y vanishes, but A takes a finite value Xa ■ 

The following exact relations between fermionic and bosonic Green's functions of original and dual variables can 
be proven: 

Goik = (A„ - Sky^g'JG^kg'J{^^ - SkY^ + (A„ - SkT^; 

(15) 

XsiK = (An - VKT'xa^nKXai^n - VkT' + (An - VkY'; 
The relation for the fermionic Green's functions is derived in Ref. IIZTl . Main idea is to substitute the right-hand side 



of Eq.( 10 1 in the formula G^jk - The last term (A^^ - ej)"' appears from the differentiation of the fermionic 

determmant. The bosonic relation is derived quite similarly, but the signs appearing during the derivation should 
be discussed. First, the expression for the bosonic Green's function has a different sign, Xoit. - \ ■ This gives 
factors of -1 at X and X. Second, since the bosonic Hubbard-Stratonovich transformation contains an inverse of 
the determinant, the same different sign appears for the last term (Aq - V^)"', so that the signs of all the terms are 
changed, compared to the fermionic relation between the real and dual fermionic Green's functions. Therefore the 
entire resulting relationship between the real and dual bosonic Green's functions obeys exactly the same form, similar 
to the results of the supersymmetric Hubbard-Stratonovich transformation |[33]| . 

It is useful to introduce the self-energy (for fermions) and polarization-operator (for bosons) corrections to the 
effective-medium approximation, defined as 

S'^G'-^'; n'=x'-^^ (16) 

We remind that Q, X are given by the expressions Q. 

The self-energy and polarization operator of the dual ensemble are given by the following expressions: 

£ = G-i-^-i; n = l-i-^-'. (17) 



Being re-expressed in these quantities, the exact relations ( 19 1 have a particularly simple form: 

(Kky'-Vk+8^ 

6 
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which resuhs in the following expression for the fermionic and bosonic Green's functions: 

(19) 

^QK = Ofn +XnnnKXnT^ + An - Va: 

One can immediately see from these formulas that the effective-medium expressions Q correspond to the Gaussian 
approximation for the dual variables: £, fl (and consequently Z', 11') vanish. 

Following the effective medium paradigm we introduce a generalized self-consistent condition for A^^ and An 
similar to the Eq.(|6]l, but for "exact" fermionic and bosonic Green's functions: 



^YjG^k, (20) 



/Tn = ^ XnK 

K 

5. Basic diagrams and simple approximations 

The simplest non-local physics can already be described in the EDMFT. In our formalism this corresponds to a 
zero-th order approximation ("no diagrams") £ = and fl = 0, and condition ^ is satisfied. We intend to preserve the 
advantages of the EDMFT, but also include the formation of collective modes and their interaction with the electrons. 
For a "purely fermionic" Hubbard-like theory this means an account of the ladder series in dual space 1281. 

For the present theory with the dual action of Eq.([TT]i, there are two basic Green's function lines (Fig. |A.l| i: a 
fermionic (full-line with arrow) and a bosonic (wavy-line) as well as two interaction vertices: a triangle (corresponding 
to the dual boson-fermion interaction - X) and a square (describing the dual-fermion interaction - y). In this paper, we 



restrict ourselves to a polarization operator, given by the electron ladder ended by two triangles (Fig. A. 2 ). Equally. 



one can consider this diagram as just a two-lines bubble, with a renormalized vertex (Fig. A. 3 ) 



fln/f - -in"^ ■ Xnjf ■ 

(21) 

Here a symbolic notation is used, so that yn is a tensor, and An is a vector with components corresponding to different 
cj, and the dots in the first line denote a scalar product vector-tensor-vector. The bare dual polarization operator (an 
empty loop of the dual Green's function) is a tensor having only diagonal components. Later we also use bare 
polarization operators for the initial variables and for the impurity model, so that 



X~0,tJ ^ f~i 

£IK — ~ hk i^uki^ai+nk+K, 



(22) 



„0,w _ 

^n - -g(jjg(jj+Q.- 

Scalar polarizabilities, respectively denoted by -^nA-' ^%.k' -^n' found by summating over fermionic fre- 

quencies, so that, for example, /V^^ = - Tjka,@ojk&u,+ak+K- 

The simplest approximation one can make is just to neglect the vertex part y. In this case /i^*^ - A^n - x^- 
Simple manipulations with the second lines of Eqs. (|4]), ( [T6| ) and ( [T8| then give 

XaK = . (23) 

A similar non-local contribution to the fermionic self-energy can be introduced to describe the renormalization of the 
fermionic Green's function G. The corresponding diagram contains one fermionic and one bosonic line (Fig. |A.5| l. 
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We will however concentrate on the analysis of bosonic quantities and for simplicity use the EDMFT result instead 
of the renormalized G propagator. This allows for a further simplification. Indeed, in the local approximation § and 
gduai jjjjfgj. jjjgjj. iQ^^i p^j.j Qjjiy ||26l[27l : 

&ajk = &a>k - gu- (24) 

Further, since we suppose that y - Q and since < p >= 0, the bosonic Green's function for the impurity problem 
becomes equal Substituting the expression (24 1 into Eq. (23 i, we obtain 



(25) 



For A = 0, this expression is clearly related to the Random Phase Approximation (RPA). To get the "textbook" 
RPA, one should completely neglect the interaction terms of the impurity problem, so that g^, = {icj - A^y^ and 
0ai,k = (ii^ - fit) The scheme with an exactly calculated can be called RPA+DMFT. In this case the RPA is 
a simple case of Fig. A. 3 To go beyond the RPA, one should be careful. As we will discuss in Sec. [6] in general 
it is wrong to work with the renormalized Green's functions but to neglect the vertex y. Because this may lead to a 
violation of the conservation laws. 

with A''^ approximated by A (that is, without putting A to 



A theory using the single-bubble approximation (21 



') is essentially the GW+EDMFT approach fT4l|. There is only a slight difference in the self-consistency condition, 
as we explain in Appendix. It is also important that for the diagram calculation of the GW-i-EDMFT approach that the 



local part of the bosonic Green's function is excluded ad hoc. As one can see from Eq. ( A. 10 1, for the EDMFT choice 
of the Bose variables, it actually does not vanish exactly. So, our consideration proposes a certain a modification 
of GWh-EDMFT. Anyhow, such modification should not strongly affect the results of the GW-i-EDMFT. The real 
advantage of our consideration is the possibility to go beyond the domain of applicability of the GW-nEDMFT, i.e. in 
is the situation when the single-bubble approximation becomes incorrect. It is considered in the next sections. 

The diagrams shown in Figs. A.4 A. 3 A. 5 form the minimal consistent set beyond the RPA to describe the 



formation of collective modes with electron-boson coupling taken into account explicitly. 



6. Charge conservation 

So far, we did not consider whether the constructed approximations conserve particle number, total spin, etc. 
However, the conservation laws are particularly important for bosonic excitations. For example, the charge (particle 
number) conservation resulting in the equation Q? < pp >sik- < jj >nK immediately leads to the requirement 
< pp >K=i)- at any finite frequency (provided that the current-current correlator < jj >n^K=o is finite). For 3D 
systems with Coulomb interaction the long-wavelength asymptotic behaviour for the density-density correlator should 
be < pp >aj^K^()'^ 7v^7)2 ' where Q„ is the plasma frequency. The RPA for free electrons is proven to obey this property, 
since ^q^^q vanishes with = {ioj - £<:) ', but this is not true for the renormalized Green's function ||36l . Indeed, 
following the standard proof |[T8l[T9] one writes A'^^ = - Yjkoj&ojk&oj+nk+K as 

^0 y/ 1 1 \ 1 (26) 

^ \lLi> - Ek - I.a> 1(C0 H- Q) - Ek+K ' '^cu+Q / + 2a; - 2„+n 



For free electrons of dispersionless E (that is, in static mean-field) the summation over a> produces a difference of 
occupation numbers, rtk - iik+K, in the nominator, that vanishes ?A K - Q. However this reasoning clearly breaks 
down for the frequency-dependent Z^, in DMFT approximation. It was known for a long time that the neglection 
of vertex corrections together with the use of renormalized Green's functions, violate charge and spin conservations 
ll36l [37l l38ll . We found that with the DMFT Green's function does not vanish at small K and therefore it follows 



form Eq. ( 25 i that EDMFT+GW is not a valid description of plasmons. All the more, plasmons do not appear in 
simple EDMFT, as one can easily observe from the second line of Eq. (jij): the local quantity + Aq cannot behave 
like Q?K^. This is the reason why the GW approach is normally used only for calculating the single-particle Green's 
function, while any two-particle properties are investigated within the Bethe-Salpeter equations [391 . Nevertheless, 
such a scheme cannot describe self-consistent effects of collective excitations on the single-particle spectral function. 
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To construct a conserving approximation, it is useful to consider the following reformulation of the theory. Let us 
return to the initial lattice action ([TJ and do the Hubbard-Stratonovich decoupling for the bosonic variables only. It 
will give an analog of Eq. ( fTO] ) with the action 

S = S^[c,c^] + p*abri + ri*abP - r]*<^b(^ - Vy^UhT], ^^7) 

Physically, S ^ describes a lattice with the interaction being local in space but retarded. Now we formally integrate over 
the fermionic variables. The Gaussian approximation in dual bosons for such a theory corresponds to the expression 

XnK = ; , (28) 

where X'^ is an exact /T-dependent two-particle correlator for the lattice problem 5^. Formally, the calculation of 
corresponds to a summation of all parts of the diagrams with fermionic hnes in series for Eq.( 111. Comparing Eqs. 
( [28] ) and ([T9| we obtain 

t^nK-Xa'x'nKXa'-Xn'- (29) 

This way, we need a conserving approximation for X'^ . Such an approximation is actually known within the 
framework of DMFT. Because it is a conserving theory in the Baym-Kadanoff sense [36 J, with the self energy being 
just the sum of all local diagrams. In such a theory, the conserving approximation for two-particle quantities can be 
obtained as a Green's function variational derivative with respect to the dispersion law, 

< (cic|, +gn')(c2cl, +822') (30) 

SS22' 

An important property of the Baym's criterion of consistency is that the variation should be considered self-consistently, 
so that the A and g of DMFT vary with e to preserve the self-consistency condition (j6|l. One can also explicitly see 
that such a self-consistent variation preserves the charge conservation. Indeed, the charge conservation follows from 
gauge invariance under the transformation Cj,- — > Cj^yC^'^'^" , cj,- — > c|,.e'*'*'' with an arbitrary dependence of the 
infinitesimal phase shift dip on time and space arguments. In an exact theory, such a transformation just results in 
similar phase rotations for averages, for example G\\i — > Gu'e'^'^'^^K A conservative approximation must inherit this 
property. DMFT is entirely defined by the self-consistency condition (j6]l and impurity action Q. If all the quantities 
entering the Green's function ff^i^ (such as s, A, g) evolve simultaneously under the gauge transformation, the self- 
consistency (|6]l clearly stays fulfilled. The entire Green's function ^^i^ obeys the gauge phase rotation, which proves 
that the theory is conserving. 

We can prove that the ladder approximation for dual bosons results in such a conserving theory. We introduce a 
variation of the local part of the dispersion law ds,,',,- - e'*^''""'*fen/i-(f - f'), and corresponding variations of local 
quantities 6A,t>r - e'^'^''^^'^6AnKit - t'), and Sg„'r = e^^^'^^^'^SgsiKit — t'). In the tensorial notation used previously, the 
variations 6snK, ^Aqk, dgnK are vectors with a single index to ; their components are defined as Fourier-transforms; 

The straightforward variation of the self-consistency condition with the DMFT Green's function (|6|) gives: 



+ 6Aqk - SsaK 



(31) 



It is worth to remind that X^^ and are tensors with double indices (w, lo'), as well as the quantities X^^ and xq to 
appear in the following formulas. 



Now, let us consider the full lattice susceptibility X^ „ as a Fourier-transform of Eq.( 30 1 with pairwise-equal space 



arguments: 
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The self-consistency requires 6Gi 



- Sgt,,t',ri, and therefore the above-introduced variations correspond to 



- - 



Similarly, for the susceptibility of the impurity model 



Sgn,K 
5Ea,K 



xn = - 



SAnK 



(33) 



(34) 



Although the right-hand side formally includes the wavevector K, x is of course /iT-independent. 
With the above two formulas formulas, it immediately follows from ( [3T| l that 



1 



(35) 



Such a self-consistent condition ( 35 1 is equivalent to locality of the irreducibly vertex and it is used in the DMFT 



scheme for calculations of non-local susceptibilities |6|. It is worthwhile to remind that the susceptibilities with a sin- 
gle frequency argument, which we defined previously, can be obtained by summation over the fermionic frequencies. 



so that X^j^ 



J oil, HIT 



and^n = 2, 



To establish a connection with the dual formalism, we substitute this expression into Eq. ( 29 1 and express X using 



Eq. (35 I. Simple algebraic manipulations give 



Xn 1 



x" - x" 
-^n^: -^n 



1 Xn 



^nxo 1 i(x„-xO)i,(x«^-x^)x0^nJ 



(36) 



One can recognize that —-5- and -ir(xn - x'')-;;- are, respectively, the previously introduced quantities A and y with 
proper indices. Furthermore, it follows from the self-consistency condition d6]) that 



"■n/r 



X° 



(37) 



So we conclude that the obtained expression for IT is exactly the result of the dual ladder summation (21 
Finally, let us consider the plasmon dispersion given by the equation 



1 



V^-An 



(38) 



In the long-wavelength limit, one considers an expansion in K\ in particular X^^ = x[^''°' + X^'^'/T^ + .... The zeroth- 
order term in the left-hand side vanishes due to charge conservation. The second-order term results in the dispersion 
equation for plasmons in a 3D correlated lattice: 



y(0,2) ,eif\ ,^ _ 



1 



(39) 



For the 2D case, the right-hand side equals 

While deriving the last equation, we supposed that An stays finite and can be therefore neglected in comparison 
with Vk at small K. Note however that this quantity still enters and therefore affects the left-hand side of the 
dispersion relation, so that the physics of local plasmonic fluctuations is taken into account. 



7. Ladder summation in a strong-coupling limit 

The single-bubble approximation becomes invalid when vertex y is large. There is a general situation of this kind: 
the four-point vertex of an impurity problem with magnetic moment formed appears to have a part displaying a slow 
dynamics that give a large contribution at low frequencies. This simply means that for isolated atoms the magnetic 
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moments are free and, thus, for almost isolated atoms (when the energy of interatomic effective interactions is much 
smaller than the intraatomic one) they can be considered as almost integrals of motion. 

To illustrate this statement, let us consider the atomic limit of the single-band Hubbard model. The single-electron 
Green's function of the Hubbard atom at half-filling 



a? + (t//2)2 



(40) 



has a time-scale tu oc which we call fast. Contrary, the spin-spin correlator < StSt' > is independent, as one can 
easily check, of time arguments. So, a time-scale for < SjSr' > is the inverse temperature /3. It means that, apart from 
a small time domain |t - t'| s; f/"', the value of < St-Sj' > is determined by its non-Gaussian part, which does not 
fall as |t - r'l increases. The much slower dynamics of y in comparison with ^ just reflects the fact that a rotation of 
the spin does not change the energy of the atom, whereas the change of the particle number does. Direct calculation 
of -y for the Hubbard atom does support this observation. The finite-temperature expressions 140) for 71234 contain a 
"singular" term proportional to pU^dno- So, the first two time arguments of y in time-domain must coincide, as well 
as the last two: y = Yttt't', but the difference r - t' can be large. 

The properties of the impurity problem are different from those of the isolated atom. However, for the strong- 
coupling limit the hybridization is small. In this case the fast dynamics roughly stay unchanged, so that the atomic 
Green's function ( 40 1 can be used, and the four-point vertex obeys the form Yttt't' ■ The hybridization can introduce 
its own slow timescale ta, so that y does depend only on t - r'. In the frequency representation, this means 



7 = rn, 



(41) 



and at small frequencies yn oc t^U^- Since y includes a large factor, we can keep only the irreducible term in the 
expressions for A andx- Note that this is completely opposite to the case considered in Section|5] 

The terms proportional to gg in Eq. ( 12i are small in and therefore can be neglected. Then we obtain the 

(42) 



following relation from Eq. (12 1; 



In the above equation we introduced a "bare susceptibility" of the impurity problem 



(43) 



Next, we use the relation ( 14 1 and keep only the large factor y: 



(44) 



Note that in the present approximation A, like y, only has single frequency argument. 

Let us consider the half-filled Hubbard lattice with the nearest-neighbour hopping f, at large U. The low-frequency 
behaviour of this system is essentially the dynamics of the Heisenberg model. The effective-medium description of 
spin waves in this case should result in the expression 



1 



+ Ao-Jr 



(45) 



with Jij - jj for nearest neighbors. Although this expression is rather simple (in our notation, it is just 11' - Jk), it 
was not reproduced by any DMFT-like approximations so far The superexchange interaction containing t^ can easily 
be obtained from a single loop of the fermionic lines, but the RPA and GWh-EDMFT denominators contain an inverse 
of this quantity (see Section|5]l, so the structure of the formulas is drastically different. 

In our formalism, the desired expression appears from the ladder summation, that is the calculation for the Eq. 



(21 1 and corresponding diagram shown in Fig. (A. 2 1. In this summation, we take into account the "slow" vertex (41 



only. Since y depends only on a single time argument, the calculation is quite simple. With the exact relation (18 1, we 
obtain 

'A) x-l 



11 



-1 



(46) 
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Finally, we substitute the formulas (42i, (44i foix and A. This results in the cancellation of y out of the expressions: 



n' 



(47) 



where An from Eq. (44i has a meaning of an effective local Stoner parameter. Both and obey the "fast" 
dynamics only. To cletermine the low-energy properties, in which we are interested, it is enough to calculate it 
at Q = 0, and to replace the summation over frequency by an integration. A straightforward calculation of the 



expression ( 43 1 with g - gat gives X%,^q ~ 4u''^° calculate ^q^q, we switch to real-space and use the strong-coupling 



limit. It gives the Green's function Q - {^g^^ t for the nearest neighbors and Q - Q elsewhere. Now, the integration 

gives Xn=o - ^ for the nearest-neighbors and we end up with the desired expression (45 1. This derivation can be 
considered as a generalization of the effective exchange formula for the magnetically (OTered case BTl |42; 1 to the 
paramagnetic state. 

The standard theory of magnons in itinerant-electron systems U is based on the RPA which is supposed to be 
exact at zero temperature; at finite temperatures, magnon-magnon and electron-magnon interactions should be taken 
into account |[Tl|2l[T2l. Electron-hole excitations are primary objects in such an approach, spin waves are coherent 
superposition of such excitations and, in general, are only well-defined in a restricted part of the Brillouin zone, 
due to Stoner damping. In the case of strong-coupling Hubbard model, we have an "anti-RPA" situation when local 
magnetic moments are primary objects. Therefore, completely diff'erent approaches should be used, e.g. those based 
on the Hubbard X-operators |4 , 43 1. It is very non trivial that the dual boson approach describes these two opposite 
cases within the same formalism, thus providing a practical way to interpolate between weak and strong coupling 
limits. 



8. Conclusions 

In summary, the general scheme for non-local correlation effects within the dual boson approach is presented. 
This method gives a direct means to obtain the higher-order diagrammatic contributions beyond the EDMFT scheme 
free from the double-counting contributions and it can be useful to describe the magnon-spectrum renormalization in 
strongly correlated systems. The method combines the numerically exact solution of an effective single-cite impurity 
model with frequency-dependent electron-electron interactions as well as an analytical diagrammatic expansion of 
non-local fermionic and bosonic self-energies. This approach is designed to describe of correlated systems with 
strong local fermion correlations and fluctuating non-local bosonic modes on equal footing. 

In this paper we focused on the applications to the description of bosonic degrees of freedom. Additional contribu- 
tions to the electron self energy also arise, in comparison with the dual fermion approach |26|, e.g., the diagram shown 
in Fig. A. 5 Similar diagrams are important in the cases of coexistence of disorder with electron-electron 1441 145) 



or electron-phonon B6I interactions. It is important to calculate these additional contributions within the Hubbard 
model and make a detailed comparison with dual fermion theory. This will be the subject of our further studies. 
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Appendix A. Dilferent choices of dual variables 

The formalism presented here has an important difference compare to the common introduction of the EDMFT: 
whereas we introduce dual bosons 77 by the Hubbard-Stratonovich decoupling of the quadratic form (V - A)p*p, the 
term decoupling Vp*p is decoupled in the standard EDMFT scheme lfT4ll47l . As one can see, this does not result in any 
physical consequences, as the two approaches yield the identical results for the averages 0,X over original variables. 
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However, properties of the dual bosons, of course, do depend on their definition. In particular, the difference is about 
local properties of the dual propagator: for the present formalism, its local part vanishes in the EDMFT approximation. 



whereas the "effective-medium" condition 



k 



0, 



(A.l) 



(A.2) 



holds for the decoupling (A. 3 i. Here ?C is the Green's function of the dual variables arising from the decoupling of 
Vp*p: 

Vp'p^ a(p*(p + fp)-V'^a^4>'^, (A.3) 
and;t'0 =< (p<p* >i,„p is defined for the impurity problem with decoupled bosonic hybridization, so that 



' imp 



(A.4) 



To maintain a similarity with our formalism, we introduced in the above expressions a scaling factor a, that is absent 
in the standard EDMFT scheme (14). All physical quantities do not depend on a particular value of this factor. 
Averages and;^^^ can be found from the exact relations similar to ( 19 1: 



a^X"" = VXV + V; 



a = AA'A + A. 



(A.5) 



We will now consider a general case of some approximation utilizing a decoupling of Wp*p with an arbitrary W. 
The approximation can be non-Gaussian, so that there is a correction 11' to the EDMFT polarization operator: 



X 



1 



The Green's function for the dual bosons can be found from the exact relation 



(A.6) 



(A.7) 



One can realize that there are two special cases. First, let us take W - V H- IT - A, and use a - x as in the main 
body of the paper. For this choice, using ( A.6 A.7| i one can straightforwardly show that 



^y+n'-A 



xiv + n-K) 



-,^^-x- 



x-^ + K-y -XV 

Summing over k and taking Eq.(|20]l into account we conclude that the condition is fulfilled 

0. 



(A.8) 



(A.9) 



It turns out that the dual propagator vanishes for variables obtained from the Hubbard-Stratonovich decoupling of 
the renormalized interaction, with the hybridization part excluded, if the condition pOb is fulfilled. For the EDMFT 



approximation, 11' is absent and the condition ( A. 1 1 holds for the decoupling with V H- 11' used in the present paper. 
Now, let us consider the case W = V -n 11'. It is suitable to take a - x^^ H- A in this case. We obtain 



X' 



4- A - y - H' 



(A. 10) 



Again taking a summation over k and substituting Eq.( 20 1, we obtain Yjk ^^^^ = A'A (a ' + . The second line of 
(A.5 I with a - x'^ + A gives the same answer for;^'^, therefore we obtain 



k 
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-■X4>- 



(A. 11) 
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The "effective-medium" condition holds for variables obtained from the Hubbard-Stratonovich decoupling of the 
renormalized interaction, with the hybridization part included, if condition ( 20 1 is fulfilled. For EDMFT, this corre- 
sponds to the decoupling of purely Vp'p part used in Ref. [14]. 



We stress that, strictly speaking, apart from the simple EDMFT case the conditions ( A. 1 1, ( A. 2 1 are not fulfilled 



simultaneously with (A.l i. This does not cause serious problems for the GW+EDMFT calculations, because one 
typically does self-consistent calculations of the EDMFT problem and then calculates diagram corrections for a pure 
EDMFT hybridization. On the other hand, for higher order diagrammatic approximations the difference can be essen- 
tial. 

Finally, we note that a specific choice of the decoupling procedure is physically irrelevant. As we have seen, 
EDMFT can be constructed using W-V-KorW-V, and the result for initial ensemble is the same. It can be 
checked that any other decoupling will also work; what matters is the Gaussian approximation for dual variables and 
the value of the hybridization. The latter is defined by the condition ( A. 1 1, that does not include any dual quantities. We 
use the theory based on W = V- A decoupling because it allows a simple construction of the dual diagrammatic series 
without a double counting of the local contributions. This important property is based on the fact that the EDMFT a 
dual Green's function has a very simple physical meaning: according to (A. 8 i, it is just the Green's function of initial 
variables with local part subtracted. 
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Figures captions: 

Fig.l Basic building blocks for dual-boson diagram: fermionic {Q) and bosonic {X) dual propagators as well as A 
and J vertices. 

Fig.2 The bosonic dual self-energy in the ladder approximation. A triangle represents the A vertex and a square 
represents the y vertex. 

Fig. 3 The bosonic dual self-energy with the renormalized triangle vertex. 
Fig.4 The diagrammatic equation for the renormalized triangle vertex. 

Fig.5 An example of diagram for the fermionic dual self-energy with the renormaUzed triangle vertex. 
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Figure A.l. Basic building blocks for dual-boson diagram: fermionic (0) and bosonic (X) dual propagators as well as A and y vertices. 



gure A.l. The bosonic dual self-energy in the ladder approximation. A triangle represents the A vertex and a square represents the y vertex 
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Figure A. 3. The bosonic dual self-energy with the renormalized triangle vertex. 
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Figure A. 4. The diagrammatic equation for the renormahzed triangle vertex. 
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Figure A. 5. An example of diagram for the fermionic dual self-energy with the renormalized triangle vertex. 
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